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Abstract 

We have analyzed archival Ginga data of a number of Soft X-ray Transient 
sources and modelled them in an attempt to constrain the geometry and physical 
properties of accretion. Our models include a self-consistent description of X- 
ray reprocessing, linking the properties of the reflected continuum and those of 
the iron fluorescence Ka line, and allowing for relativistic smearing of spectral 
features. We derive constraints on the inner extent of the disk and its ionization, 
assuming a generic geometry of a central X-ray source and an external accretion 
disk. 

Fitting our models to the data we have found that in some sources - 
GS 1124-68 (Nova Muscae 1991) and GS 2000+25 - the evolution during de- 
cline proceeded in qualitative agreement with recent theoretical models based on 
advection-dominated solutions of accretion flow. These models link the high/soft- 
low/hard transition with a change of geometry, most importantly, the inner radius 
of the optically thick disk. Quantitatively, our analysis requires serious revision of 
the model assumptions. A remarkably exceptional case is that of GS 2023+338, 
where the source behavior both during outburst and decline was quite different 
than the other, more "standard" ones. 

1. Introduction 

Accretion flow around black holes in compact objects (both active galaxies 
and stellar sources) clearly proceeds through at least two phases: optically thick, 
cooler plasma producing a soft thermal component, and optically thin, hot plasma 
where hard X-rays/7-rays are produced (see Mushotzky et al. 1993 and Tanaka 
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& Lewin 1995 for reviews). 

The geometrical configuration of the two phases is still uncertain. Simple 
possibilities include pure radial and pure vertical stratification of the flow, but 
more complicated geometries, with partial overlap are also considered. Localized 
active regions above the disk, appearing due to e.g. magnetic activity, are also 
considered (Galeev et al. 1979; Stern et al. 1995). 

The commonly made assumption that the observed hard X-ray spectra 
are due to the inverse Compton process requires an input of soft photons from 
the optically thick plasma to the optically thin one. The reverse interaction 
must then also take place: the hard X-rays will illuminate the optically thick 
plasma (Guilbert & Rees 1988) and they will be reflected and reprocessed. As 
a result, a continuum spectral component is formed by Compton reflection, with 
spectral features due to emission/absorption by heavy elements (primarily iron) 
superimposed on it. 

Studying the effects of such reprocessing can be extremely useful as a tool 
for constraining the spatial distribution of the two phases of accreting plasma. 
The amplitude of the reflected component constrains the solid angle subtended 
by the optically thick plasma from the X-ray source. The shape of the reflected 
continuum below ~ 10 keV, and the energy and strength of iron spectral features 
around 7 keV give information about ionization state of the reflecting plasma. 
Possible broadening/smearing of spectral features can constrain the plasma ve- 
locity field if it can be distinguished from e.g. broadening due to Comptonization. 

Black hole Soft X-ray Transients (SXT) are particularly suited to this 
kind of analysis. Those are a sub-class of X-ray binaries that undergo occasional, 
dramatic outbursts. During the subsequent decline phase the mass accretion 
rate, which is thought to be the most important parameter, decreases roughly 
exponentially, with a typical e-folding time scale of 30 — 60 days. During the 
decline they go through a sequence of well defined states characterized by distinct 
spectral and timing properties (review in Tanaka & Shibazaki 1996). 

2. Data 

For our analysis we use archival Ginga data of a number of SXT. The data 
combine broad band (1 to 20-30 keV) with moderate spectral resolution (18% at 
6 keV), very good statistics due to large effective area of the LAC detector (~ 
4000 cm~ 2 ) and systematic uncertainties of the instrument below 1%. Background 
subtraction can pose a problem for bright sources, due to the contamination of the 
background monitors by source counts. We have developed a method to estimate 
the background, described in some detail in Zycki et al. (1998b). 
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Fig. 1. X-ray reprocessed spectra showing effects of ionization (dashed vs. 
dot-dashed curves) and relativistic smearing (dotted vs. dot-dashed curves). 



3. Models 

Our model of X-ray reprocessing, described in details in Zycki et al. (1998b), 
self-consistently combines the model for Compton-reflected continuum and Monte 
Carlo computations of the iron Ka line. The main model parameter is the ion- 
ization parameter, £ = AirF/n, where F is the illuminating flux and n electron 
density. The amplitude of the component is expressed as a solid angle of the 
reflector, normalized to 2tt, f = Q/2tt. The effects of relativistic smearing are 
computed convolving the total reprocessed spectrum with relativistic line profiles 
implemented in XSPEC as the "diskline" model. The main parameter here is 
the inner disk radius in units of R g = GM/c 2 , assuming a fixed illumination 
emissivity (we assume F iTT (r) oc r~ 3 ). Figure [l] shows typical model spectra. 

4. Results 



4-1. Low/hard state 

The results are shown in Table |T] The reprocessed component is signifi- 
cantly present in the spectra. Its amplitude is significantly smaller than 1, typical 
values being 0.3 - 0.7. It is rather weakly ionized, £ < 50 corresponding to mean 
ionization level of iron <FeXVII i.e. only M-shell electrons removed. Relativistic 
smearing of the reprocessed component is significant in some of the spectra, and is 
consistent with decreasing in time, as the sources progressed through the decline 
phase. Even at the beginning of decline in GS 2023+338 the smearing is less than 
expected from a disk extending down to the last stable orbit at r ms = 6 R g , typical 
values b eing -R; n ~ 20 30 -Rg (see Figure ^] and Zycki et al. 1997), although R; n 
is compatible with r ms in Nova Muscae (Zycki et al. 1998a). 
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Table 1. Results of fitting the reprocessing model to the data 
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4-2. High/soft state 

In the high state spectra the major difference in the reprocessed compo- 
nent is its strong ionization (Table |Tj). Dominant iron ions are usually He- and 
H-like ones. Consequently, the Ka line and absorption edges are much stronger 
in the composite spectra. The amplitude / is also larger, close to 1. Relativis- 
tic smearing is significant, but the values of _R in differ between different objects. 
The fact that, for GS 2000+25, R in ~ 6 (which is not expected given the as- 
sumed -Fin- oc r~ 3 i.e. without the boundary condition term; Shakura & Sunyaev 
f 973) may suggest that an additional broadening of that component (e.g. due to 
Comptonization) takes place. 

We find that the soft component of the continuum spectrum cannot be 
described by either pure black body or multi-temperature (disk) blackbody spec- 
trum. Additional power law tail is required, which can be modelled as additional 
Comptonization of a (disk) black body. This differs from results of Ebisawa et 
al. (f 994) who were able to fit the disk blackbody to the spectra but they were 
using the smeared edge model with a narrow gaussian to model the iron spectral 
features, rather than a reprocessing model. 

4-3. Very high state 

Nova Muscae showed an example of Very High State spectrum at the peak 
of its outburst (Jan-Feb 1991). On 11th Jan the soft continuum component 
again cannot be modelled as a (disk) blackbody, so we assume a Comptonized 
blackbody instead (Figure |3], middle panel). The presence of the reprocessed 
component is then again highly statistically significant although its amplitude is 
not well constrained from above, in particular if radial distribution of ionization 
is allowed for. The reflection is strongly ionized and smeared, similarly to high 
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Fig. 2. Spectra of Nova Muscae in low and high states (upper panels) and \ 2 
contributions (lower panels). Lower solid histograms show the iron Ka line and 
the reflected continuum, dotted histograms show the soft thermal component 
and the dashed histograms show the primary power law. See Table [T] for fits 
details. 



state spectra (Table [[]). Thus although the overall spectrum is dominated by the 
hard power law (r ~ 2) as in LS, the properties of both the soft component and 
the reprocessing are similar to HS, so the suggested histeretic behavior of SXT 
(Miyamoto et al. 1995) does not seem to be complete. 

5. The exceptional case of GS 2023+338 

GS 2023+338 was clearly an unusual source, showing fast and chaotic 
flux and spectral variability, and no clear soft, thermal component even when its 
luminosity was very high (Tanaka & Lewin 1995). 

The spectrum of GS 2023+338 at the peak of its outburst was rather un- 
usual: it can be described by an optically thick comptonization of a disk blackbody 
radiation of rather high temperature, kTo fa 1.4 keV in a rather cool plasma cloud, 
kT e fa lOkeV, r fa 6. Comparison with spectra of Nova Muscae at a similar lumi- 
nosity level (Figure |3]) suggests that the actual mass accretion rate in GS 2023 
might have been super-Eddington. The hypothesis that the observed emission is 
a (quasi-thermal) disk radiation is not however easily compatible with the pres- 
ence of spectral features due to highly ionized iron, since in such a geometry disk 
irradiation is rather ineffective. 

The short time-scale (~few sec) flux and spectral variability can be at- 
tributed to strong and rapidly variable photo-electric absorption. During one the 
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Fig. 3. Comparison of two spectra of Nova Muscae in Very High State and a 
spectrum of GS 2023+338 at the very peak of its outburst. The luminosities are 
bolometric ones but the correction is rather uncertain in the case of GS 2023. 



episodes of such strong absorption on 30th May, the source spectrum is consistent 
with a presence of a soft thermal component of temperature T « 1 keV, and a 
steep power law tail, r » 2. Such spectra are typically observed in high/soft 
states of GBH (Figure |]; Zycki et al. 1998c). 

During its decline from outburst the source showed a fairly regular, expo- 
nential decline of the luminosity (to ~ 30 days) with again rapid variability due to 
photo-electric absorption superimposed on it. The amplitude of the reprocessed 
component decreases with time, in rough correlation with decreasing level of rel- 
ativistic smearing, i.e. the behavior is consistent with R in increasing. However, 
the power law spectral index remains practically constant, in sharp contrast to 
the behavior of GS 1124-68, GS 2000+25 and GX 339 (Zycki et al. 1998b). 

6. Discussion 

The results presented in previous Sections clearly suggest that the geome- 
try of accretion and the physical state of reprocessing material change as sources 
make transition between hard and soft states. The most significant changes to the 
reprocessing are those of ionization of the reflector and the amplitude of reflection. 
In the low/hard state the solid angle of the reprocessor from the X-ray source is 
smaller and the reprocessor is colder than in the high/soft state. The amount of 
relativistic smearing does not show clear, statistically significant changes between 
the transitions, but i? in is consistent with increasing as the amplitude / decreases. 

We can compare the above results with predictions of models proposed for 
the spectral transitions. 
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Fig. 4. Power density spectra in three energy bands (a) and energy spectrum (b) 
of GS 2023+338 on 30th May (around 7:00 GMT), suggesting that the source 
was in the high/soft state, but it coincided with heavy absorption. 

The observed changes are qualitatively expected in the models of Mineshige 
(1996) and Esin et al. (1997). These models assume predominantly radial stratifi- 
cation of the flow, with optically thick disk extending only outside certain transi- 
tion radius. The central, hot solution derives from advection dominated accretion 
flow solutions although when applied to typical low state, it has radiative efficiency 
of ~ 0.35. The transition radius (which can be identified with our R m ) has to de- 
crease with increasing m to explain softer spectra when sources are brighter, and 
this also explains larger amplitude of the reprocessed component and its strong 
ionization in high state. The behavior of R in we find from relativistic smearing is 
also consistent with that prediction. Thermal radiation emitted by the disk can 
by Comptonized by hot upper layer of the disk, giving rise to the observed steep 
power law tail. 

Quantitatively, however, the values of R in we derive are rather smaller 
than those postulated by Esin et al. In particular, in the low state we derive 
R in ~ 20 - 50 R g , while they assume 2 x 10 4 R g . Esin et al. (1998) argue that the 
smaller values are still compatible with their model but this claim has not been 
verified by actual data modelling. 

Relatively small values of i?i n we derive indicate that there must exist 
significant overlap between the central hot X-ray source and external cool disk 
i.e. the hot source has to extend much further out than R in : otherwise the energy 
release within the source would be much too small to account for the overall 
energetics and spectral shape (Gierlihski et al. 1997; review in Poutanen 1998). 
Perhaps the hot plasma forms a corona at large distances. The transition between 
the disk and the corona at large distances could result from e.g. the thermal 
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instability of X-ray irradiated plasma, as considered by Witt et al. (1997). 

An alternative geometry was proposed by Di Matteo et al. (1998). They 
suggest that magnetic flares above the disk produce the hard X-ray radiation, the 
disk itself always extends to the last stable orbit, and the hard/low - soft/high 
transition is due to changes in characteristic height of the flares. In order to 
account for harder spectra in the low state they postulate that the flares are 
higher above the disk in this state. This however does not naturally explain why 
the reflection amplitude is smaller in the hard state unless the flares are strongly 
concentrated towards the inner edge of the disk, so a significant fraction of photons 
is lost within the central hole. An alternative explanation, that the inner disk is 
strongly ionized and so the iron spectral features are weak, does not seem to be 
supported by data analysis, at least in Cyg X-l (Done & Zycki 1998). 
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